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We present a new class of boron sheets, composed of triangular and hexagonal motifs, that are 
more stable than structures considered to date and thus are more likely to be the precusors of 
boron nanotubes. We describe a simple and clear picture of electronic bonding in boron sheets 
and highlight the importance of three-center bonding and its competition with two-center bonding, 
which can also explain the stability of recently discovered boron fullerenes. Our findings call for 
reconsideration of the literature on boron sheets, nanotubes and clusters. 
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All boron nanotubes (BNT), regardless of diameter or 
chirality, are predicted to be metallic and have large den- 
sities of states (DOS) at their Fermi energies (E^?) [1 . In 
contrast, carbon nanotubes (CNT) can be semiconduc- 
tors or metals with small DOS at their E^. Metallic CNT 
are used widely to study one-dimensional (ID) electron- 
ics [2j [3] and are superconducting at low temperatures 
H [5]. Due to the larger DOS, BNT should be better 
metallic systems for ID electronics and may have higher 
superconducting temperatures than CNT. 

Recent experiments have fabricated boron nanotubu- 
lar structures both as small clusters [6 and long, ID 
geometries [7 . Understanding the properties of BNT is 
crucial for realizing their applications. For CNT, it has 
been fruitful to study two-dimensional (2D) graphene: 
e.^., many properties of CNT are derived directly from 
graphene [3 [9]. For boron, however, no 2D planar struc- 
ture exists in its crystals which are built from B12 icosa- 
hedra [10 . Researchers have proposed several 2D boron 
sheets (BS). The hexagonal graphitic BS was found to be 
unstable [TT1[T8]. Based on extensive theoretical studies 
of boron clusters [TTJ |T2] [131 HH US] 5 Aufbau principle 
was proposed whereby the most stable structures should 
be composed of buckled triangular motifs [12 . Experi- 
ments on small clusters of 10-15 atoms support this view 
[T6^. A recent study of many possible sheet structures 
found, again, the buckled triangular arrangement to be 
most favorable [17]. Hence, 2D triangular BS have been 
studied and used to construct BNT |18l [T9l [20] . 

In this Letter, we present a class of boron sheets that 
are more stable than the currently accepted ones. We de- 
scribe their structures, energetics, electronic states, and 
provide a clear picture of the nature of their bonding that 
clarifies their stability. We also show that clusters with 
these structures are competitive with or more favorable 
than those considered to date. Our findings have im- 
portant consequences for understanding and interpreting 
the properties of these systems. For example, the un- 
usual stability of Bgo fullerenes [21] can be explained by 
our bonding picture. Hence, in our view, it is necessary 
to reconsider previous work in this general field. 

We use Density Functional Theory [22l [23] within the 
ah initio supercell planewave pseudopotential total en- 



ergy approach [2T. Calculations are done by PARATEC 
[25^ . We use both the local density approximation 
(LDA) [23l [26] and the generalized gradient approxima- 
tion (GGA) [2T for exchange and correlation. Most re- 
sults below employ the LDA and key results are checked 
by the GGA. The LDA and GGA yield same qualitative 
results with minor quantitative variances. The planewave 
basis has a 32 Ryd cutoff energy. K-point samplings 
for each case converge total energies to better than 1 
meV/atom. Norm-conserving pseudopotentials have cut- 
off radii r^=1.7 and r^=2.1 a.u.. The BS are extended in 
x-y directions while supercells have periodic copies along 
z where a separation of 7.4 A is sufficient for conver- 
gence. For all structures, relaxations are preformed un- 
til the atomic Hellmann-Feynman forces are less than 1 
meV/A and all in-plane stresses are less than 5 MPa. 

Table [Tj shows our results for four different sheets: the 
flat and buckled triangular sheets [181 Il2, , the hexagonal 
sheet, and one of our sheets (a in Figure [T]). The hexag- 
onal sheet is unstable with respect to in-plane shear, so 
we obtain the values in Table [l] by maintaining hexag- 
onal symmetry while optimizing the bond length. The 
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FIG. 1: (A, B) Two examples of our BS (top view). Red solid 
lines show the unit cells. (C) Four boron clusters: B24(a) and 
B32(a) are clusters with hexagonal holes; B24(^) and 632(6) 
are the double-ring clusters from refs. [T4l[T5]. Gray balls are 
boron atoms, and gray "bonds" are drawn between nearest 
neighbors. 



2 



TABLE I: Binding energies Eb (in eV/atom) and geometric 
parameters (in A) of four BS: the flat and buckled triangular 
sheets, the hexagonal sheet, and one of our sheets (a in Figure 
[l]). d^^"^^ is the bond length of the flat triangular sheet, d"^ 
and G?^*"^ are the bond lengths of the buckled triangular sheet, 
d^ is between atoms with the same z, while c?^*"^ is between 
atoms with different z. Az is the buckling height, d^^^ is 
the bond length for the hexagonal sheet, d^^^ gives the bond 
length range in the sheet a. 
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"Boron's atomic spin-polarization energy of 0.26 eV/atom ex- 
plains the Et differences between [9] and our work or |18| . 

^While the absolute Et from FTTI do not match our GGA results, 
Et differences among the sheets match very well. 



binding energy is 

Eb = Eat ~ E sheets 

where Eat is the energy of an isolated spin-polarized 
boron atom and E sheet is the energy per atom of a sheet. 
The buckled triangular sheet is more stable than the flat 
one due to the former forming stronger cr bonds along the 
buckled direction [T9^. We also can reproduce previous 
results on BNT made from triangular sheets [18l [19] . 

Figure [T] shows two examples of our BS, which are more 
stable than the buckled triangular sheet by 0.12 (a) and 
0.08 {(5) eV/atom. All our sheets are metallic, flat and 
composed of mixtures of hexagons and triangles. Sheet 
a is the most stable structure in our library. 

Our sheets can be obtained by removing certain atoms 
from a flat triangular sheet. Each removal produces a 
hexagonal hole, generating a mixture of hexagons and 
triangles. We define a "hexagon hole density" 

No. of hexagon holes 
No. of atoms in the original triangular sheet 

The triangular sheet has 7^=0, the hexagonal 77=1/8, and 
sheets a and have r] of 1/9 and 1/7, respectively. 

A priori^ the energies of these sheets can depend on 
both 1] and the pattern of hexagons. This results in a 
huge phase space of hexagonal patterns for a given r]. 
The most stable structures occur when the hexagons are 
distributed as evenly as possible at fixed r]. Figure |2] 
shows the LDA binding energies E^ versus rj for this class 
of structures. E^ reaches a maximum of 6.86 eV/atom 
at r]=l/9 (sheet a). We also have investigated the other 



extreme where hexagons form lines {e.g., sheet p). These 
"linear" structures are more stable than the buckled tri- 
angular sheet for r] ~l/9 but are less stable than the 
"evenly-distributed" class described above. 

To explain the stability of these structures, we de- 
scribe the nature of their electronic bonding. Generally, 
in-plane bonds formed from overlapping sp"^ hybrids are 
stronger than out-of-plane 7r-bonds derived from Pz or- 
bitals, so a structure that optimally fills in-plane bonding 
states should be most preferable. Guided by this princi- 
ple. Figure [3] shows projected densities of states (PDOS) 
for five BS with separate in-plane (the sum of Px and 
Py) and out-of-plane {pz) projections. 

We begin with the hexagonal sheet, a textbook sp"^ 
bonded system. All sp"^ hybrids are oriented along near- 
est neighbor vectors so that overlapping hybrids produce 
canonical two-center bonds. A large splitting ensues be- 
tween in-plane bonding and anti-bonding states. The Pz 
orbitals form their own manifold of bonding and anti- 
bonding states. The Pz PDOS vanishes at the transition 
point between the two. In the case of graphene, the four 
valence electrons per atom completely fill the sp"^ and the 
Pz bonding states, leading to a highly stable structure. 
However, a boron atom has only three valence electrons. 
As shown in Figure [sj some of the strong in-plane sp'^ 
bonding states are unoccupied, explaining the instabil- 
ity of this sheet. For our discussion below, this sheet is 
highly prone to accepting electrons to increase its stabil- 
ity should they be available from another source. 

Next, we consider the flat triangular sheet. Each atom 
has six nearest neighbors but only three valence elec- 
trons. No two-center bonding scheme leads to a proper 
description. Previous work has noted qualitatively that 
a three-center bonding scheme exists here [T^. We now 
present a detailed model of the three-center bonding with 
crucial implications for the stability of our sheets. Fig- 
ure [4] shows a choice of orientations for the sp'^ hybrids 
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FIG. 2: LDA Et v.s. hexagon hole density 77 for sheets with 
evenly distributed hexagons. The dashed line shows Et for the 
buckled triangular sheet. The solid blue curve is a polynomial 
fit. The two limiting cases 77 = and 77 = 1/3 correspond to 
the flat triangular and hexagonal sheets, respectively. Maxi- 
mum Et occurs for sheet a (77 = 1/9). 
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where three hybrids overlap within an equilateral trian- 
gle formed by three neighboring atoms. For an isolated 
such triangle, we have a simple 3x3 tight-binding prob- 
lem with symmetry. Its eigenstates are dictated by 
group theory: one low-energy symmetric bonding orbital 
h and two degenerate high-energy anti-bonding orbitals 
a*. (This is "closed" three-center bonding; details on this 
and other types of bonds are found in standard references 
[28].) These orbitals then broaden into bands due to 
inter-triangle couplings. Separately, the Pz orbitals also 
broaden into a single band (not shown). In Figure |3j the 
in-plane PDOS becomes zero at the energy separating in- 
plane bonding and anti-bonding states. Ideally this sheet 
would be most stable if: (i) two electrons per atom would 




E(eV) 

FIG. 3: Projected Densities of States (PDOS) for four BS. 
DOS are projected to in-plane (sum of s, px and py) and out- 
of-plane orbitals {pz)- Red solid lines show in-plane and blue 
dashed lines show out-of-plane projections; thick vertical solid 
lines show E f ■ AH curves are broadened using Gaussians with 
a width of 0.3 eV. The vertical scale is arbitrary. 



completely fill the 6-derived in-plane bonding bands, (ii) 
the anti-bonding a* -derived bands were empty, and (iii) 
the remaining electron per atom would half fill the low- 
energy (bonding) portion of the p^-derived band. This 
would mean that the Ei? would be at the zero point of 
the in-plane PDOS in Figure [3) Clearly, this picture is a 
valid zeroth-order description. However, E^? lies slightly 
above the ideal position and makes some electrons occupy 
in-plane anti-bonding states. In other words, this sheet 
prefers to donate these high-energy electrons which has 
critical implications below. (Although we seem to break 
symmetry by making half of the triangles filled and half 
empty, filling the entire 6-derived in-plane bonding band 
makes all hybrids equally occupied. This restores full in- 
plane symmetry: i.e., the two possible initial orientations 
of hybrids give the same final state.) 

The fiat triangular sheet, however, buckles under small 
perturbations along z [18 . The buckling mixes in-plane 
and out-of-plane states and can be thought of as a sym- 
metry reducing distortion that enhances binding. As 
shown in Figure [3} some states move below the Ei? as 
indicated by the small peak immediately below the E^?. 

Finally, we turn to the new structures. The above 
discussion has shown that the hexagonal sheet should 
be able to lower its energy by accepting electrons while 
the fiat triangular structure has a surplus of electrons 
in anti-bonding states. From a doping perspective, the 
three-center fiat triangular regions should act as donors 
while the two-center hexagonal regions act as acceptors. 
Thus if the system is able to turn into a mixture of 
these two phases in the right proportion, it should be 
able to benefit from the added stability of both subsys- 
tems. Specifically, the hexagon-triangle mixture with the 
highest stability should be the one that places the Ei? 
precisely at the zero-point of in-plane PDOS, filling all 
available in-plane bonding states and none of the anti- 
bonding ones. The remaining electrons will fill the low- 
energy p^-derived states, leading to a metallic system. 
These expectations are born out clearly in Figure |3] as 
well as in the energetic stability of the various structures 
(Figure [2|. In fact, the most stable sheet, a, satisfies this 
condition precisely while the less stable sheet, has a 
slight shift of Ei? from the ideal position. 

These findings have ramifications for boron clusters. 
Our structures and bonding picture can explain that the 
extreme stability of Bgo fullerenes composed of triangu- 
lar motifs with pentagonal holes [21] is due to the well 




FIG. 4: Three-center bonding scheme in flat triangular 
sheets. Left: orientation of sp^ hybrids. Center and right: 
overlapping hybrids within a triangle (D3 symmetry) yield 
one bonding (6) and two anti-bonding (a*) orbitals. These 
then broaden into bands due to inter-triangle interactions. 
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balance of three-center and two-center bonds. Also the a 
sheet can be seen as the precursor of Bgo just as graphene 
is the precursor of carbon fullerenes. We also have stud- 
ied some clusters. Figure [l] shows the double-ring struc- 
tures for B24 and B32 [HI [15 along with clusters con- 
structed by us. The new B24 cluster with a hexgon hole 
is less favorable by 0.08 eV/atom while the B32 is more fa- 
vorable by 0.03 eV/atom than the corresponding double- 
ring. The stability of our sheets, of Bgo, and our clus- 
ters with hexagonal holes suggests that for boron systems 
with more than 20-30 atoms, the Aufbau principle breaks 



down and a more general structural rule is required. 

In brief, we demonstrate a novel bonding mechanism in 
pure boron compounds arising from the competition be- 
tween two- and three-center bonding. This explains the 
stability of our boron sheets as well as larger boron clus- 
ters. Our results have relevant implications on the sta- 
bility and structure of boron clusters, boron nanotubes, 
and other boron systems. 

We thank the Bulldog parallel clusters of the Yale High 
Performance Computing center for providing computa- 
tional resources. 
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